BIHEP-th-01-01 
KOBE-FHD-01-02 
FUT-01-01 
August 2001 

Extraction of Polarized Gluon Distributions 
from Large— Light Hadron Pair Production 

Yu-Bing Dong|] 

Institute of High Energy Physics, Academia Sinica, 
Bezjmg 100039, P. R. Chma 

T. Moriif] 

Faculty of Human Development and 
Graduate School of Science and Technology, 
Kobe University, Nada, Kobe 657-8501, Japan 

and 

T. YamanishiQ 

Department of Management Science, 
Fukui University of Technology, Gakuen, Fukui 910-8505, Japan 



Abstract 

We propose a new formula for extracting the polarized gluon distribution from the large-p^ 
light hadron pair production for semi-inclusive processes in polarized deep inelastic scattering. In 
general, a large-pr hadron pair is produced via photon-gluon fusion(PGF) and QCD Compton 
at the lowest order of QCD. The PGF gives us a direct information on Ag in the nucleon, 
while QCD Compton becomes background to the signal process for extracting Ag. We show 
that the contribution from this background, i.e. the QCD Compton process, can be removed by 
using symmetry relation among fragmentation functions and taking an appropriate combination 
of various light hadron pair production processes, and thus the double spin asymmetry can be 
described in terms of Ag/ g alone. 
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1 Introduction 



The nucleoli is not an elementary but compound particle. Accordingly, its spin is carried by its 
constituents and is described by a sum rule, 

l = iAS + A(7 + L, , (1) 

where | on the left hand side denotes the spin of a longitudinally polarized nucleon, while AS, 
A(yf and L~, on the right hand side represent the amount of the spin carried by the quarks, gluons 
and their orbital angular momenta, respectively. In the static quark model with SU(6) symmetry, 
the nucleon's spin is totally carried by quarks alone and thus, each term on the right hand side 
of eq.(|l|) becomes AS = 1, A^f = and = 0. Moreover, in the naive quark-parton model, the 
value of AS is predicted to be ~ 0.7, which means that the quarks almost carry the nucleon's 
spin||l]|. However, triggered by the first measurement by the EMC in 1988 of the longitudinally 
polarized structure function of the proton g\ for wide kinematic range of Bjorken x and Q^, the 
situation has drastically changed, that is, the quarks carry the proton spin little and the strange 
quark is negatively polarized to the proton spin[0. These results are in significant contradiction to 
the traditional understandings based on the static quark model and/or the quark-parton model. 
Since then, a great deal of efforts have been made experimentally and theoretically for disclosing 
the origin of the nucleon spin[Q. So far, based on the next-to-leading order QCD analyses on the 
X and dependence of the longitudinally polarized structure function gi, we have got a number 
of excellent parametrization models of the polarized parton distribution functions (pol-PDFs) 
from fitting to experimetal data, with high precision, of polarized deep inelastic scattering (pol- 
DIS) for various targets[P-0. All of these parametrization models tell us that AS is around 
0.3 or smaller, that is, the contribution of the quark spin to the nucleon spin is rather small. 
The remainder should be compensated by and L^. To understand the physical ground of 
these results is the main subject of the so-called proton spin puzzle, which is still challenging 
problem to be solved. To solve this problem, it is very important to precisely know how the 
gluon polarizes in the nucleon. It is well-known that the next -to-leading order QCD analyses 
on gi bring about information on the first moment of the polarized gluon A(7|^. However, there 
are large uncertainties in Ag extracted from gi alone. Knowledge of /S.g is still limited because 
it is very difficult to directly obtain such information from existing data. 

So far, a number of interesting proposals such as direct prompt photon production in polarized 
proton-polarized proton collisions 0, open charm|]lO| and Productions in polarized lepton 



scattering off polarized nucleon targets, were presented for studying longitudinally polarized 
gluon distributions. Recently, HERMES group at DESY reported the first measurement of the 
polarized gluon distribution from di-jet analysis of semi-inclusive processes in pol-DlS, though 
only one data point was given as a function of Bjorken x[r2|. In general, a large-pr hadron pair 
is produced via photon-gluon fusion (PGF) and QCD Compton at the lowest order of QCD (Fig. 
1). The PGF gives us a direct information on the polarized gluon distribution in the nucleon, 
while QCD Compton becomes background to the signal process for extracting the polarized gluon 
distribution |ll3|, When we consider only the case of heavy hadron pair productions, we could 
safely neglect the contribution of QCD Compton processesfl^ because the contents of heavy 
quarks in the proton are extremely small and furthermore the probability of fragmenting of light 
quarks to heavy hadrons is also small. However, in this case the cross section itself is rahter small 
at the energy in the running(HERMES) or forthcoming(COMPASS) experiments and thus, we 
could not have enough data for doing detailed analysis. Furthermore, if we consider only small 
X region where the PGF dominates over the QCD Compton, we could also extract the polarized 
gluon distribution without difficulty [Q. However, here we are interested in the polarized gluon 
distribution for wide reigon of x to calculate the first moment of the polarized gluon. For light 
hadron pair productions in wide x region, QCD Compton contribution is not necessarily small 
and hence, it is rather difficult to unambiguously extract the behavior of the polarized gluon 
from those processes. 

In this work, getting over these obstacles, we propose a new formula for clearly extracting 
the polarized gluon distribution from the light hadron pair production of pol-DIS by removing 
the QCD Compton component from the cross section. As is well-known, the cross section of 
the hadron pair production being semi-inclusive process, can be calculated based on the parton 
model with various fragmentation functions. Then, by using symmetry relations among frag- 
mentation functions and taking an appropriate combination of various hadron pair production 
processes, we can possibly remove the contribution of QCD Compton components from the cross 
section and thus, get clear information of the polarized gluon distribution from the remaining 
PGF components. Here, to show how to do this practically, we consider the light hadron pair 
production with large transverse momentum as shown in Fig. 1. 



2 Cross sections and double spin asymmetry for large— 
pion pair production 

Let us consider the process of i + N ^ i' + hi + h2 + X in polarized lepton scattering off polarized 
nucleon targets (Fig. 1), where hi and /i2 denote light hadrons in a pair. As mentioned above, 
the spin-dependent differential cross section at the leading order of QCD can be given by the 
sum of the PGF process and QCD Compton as follows, 

dAa = dAapGF + dAaqco ■ (2) 

with 

dAa = da^+ — da^^ + da^_ — da++ . 

Here, for example, da-+ denote that the helicity of an initial lepton and the one of a target 
proton is negative and positive, respectively. Each term on the right hand side of eq.(^ is given 

by 

dAapGF ~ Ag{r],Q')dAapGF E e^j^f^^i, g')I^^(^2, Q') + (1 2)} , (3) 

i=u,d,s,u,d,s 

dAaQCD ~ E _ e?A/,(r/,Q2)rfAaQcz){/^J^(^i,g')/^;^(^2,g') + (l-2)} , (4) 

q=u,d,s,u,d,s 

where Ag{ri,Q^), Afq{ri,Q^) and D^i^z^Q"^) denote the polarized gluon and g-th quark distri- 
bution functions with momentum fraction rj and the fragmentation function of a hadron h with 
momentum fraction z emitted from a parton z, respectively. dAapcr and dAaqcn are the po- 
larized differential cross sections of hard scattering subprocesses for ^g — > t'qq and ^ 'q^^ i' 5'V 
at the leading order QCD, respectively. 

Here we consider the following 4 pairs of a combination for the produced hadrons hi with zi 
and /i2 with Z2, 

(i) (7r+,7r-), (ii) (7r-,7r+), (iii) (7r+,7r+), (iv) 

where (particle 1, particle 2) corresponds to {hi with zi, /i2 with Z2). Then, the differential cross 
section of eq.(0) for each pair can be written as 

(i) 

rfAa'^^'^" ~ Ag{r^,Q')dAapGF{lDf{zi,Q')Df{z2,Q') + lDf{zi,Q')Df{z2,Q') 
+ l-Df{zi,Q')Df{z2,Q') + {7r^{zi)^7r-{z2)) } 



+ ^Auiv, Q')dAaQCD {Df {z,, Q')Df {z2, Q') + {z,, Q')Df {z,, Q')} 
+ (contributions from Ad, As, Au, Ad and As) , (5) 



rfAa'^--^ ~ Ag{r^,Q')dAapGF{lD:{z,,Q')Df{z2,Q') + lDl{z,,Q')Df{z2,Q') 

+ \D:-{z,,Q')Dt{z,,Q') + {7:-{z,)^Ti+{z,)) } 

+ ^Au{^, Q')dAaQCD {Df {z„ Q')Df {z„ Q') + of {z^, Q')Df {z^, Q')} 

+ (contributions from Ad, As, Au, Ad and As) , (6) 



m 



dAa--^-^ ~ Ag{r^,Q')dAapGF{lDf{z„Q')Df{z„Q') + ^Df{z„Q')Df{z2,Q') 

+ ^-Df{z„Q')Df{z„Q') + (vr+(^0 ^ 7r+(.,)) } 

+ ^Au{r^, Q')dAaQCD [di\z^, Q')Df {z^, Q') + Df{z2, Q')Df [z^, Q')} 

+ (contributions from Ad, As, Au, Ad and As) , (7) 

(iv) 

dAa---- ~ Ag{i^,Q^)dAapGF{\Dl{z,,Q^)Dl{z,,Q^)+^-Dl{z,,Q^)Df{z2,Q^) 

+ ^-D:~'{z,,Q^)DI-{z2,Q^) + (7r-(zi) ^ 7r-(z2)) } 

+ ^Au{r^, Q')dAaQCD {Df (z, , Q')Df (z^, Q') + (z^, Q')Df (z, ,Q')} 

+ (contributions from Ad, As, Au, Ad and As) . (8) 

Based on the isospin symmetry and charge conjugation invariance of the fragmentation functions, 



various fragmentation functions in eqs.(^)-(P) can be classified into the following 4 functions[16|, 
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where D and D are called favored and unfavored fragmentation functions, respectively. Consid- 
ering the suppression of the s quark contribution to the pion production compared with the u 
and d quark contribution, we do not identfy Dg with D. This seems to be confirmed by 'leading 
paritcle' measurements [|l^, |T8|. By using these 4 kinds of pion fragmentation functions, we can 



make an interesting combination of cross sections which contains only the PGF contribution as 
follows; 

rfAa"^"" + rfAfx"""^ - t/Afx""'"^ - rfAa"""" ~ ^^^^^^ Q^)dAapGF (9) 



X 



From this combination, we can calculate the double spin asymmetry All defined by 



At T — 

Ag{7],Q^) dAapGF ^^^^ 



giv, Q"^) dapGF 

where the factor of the fragmentation function in eq. (|^) is dropped out from the numerator and 
the denominator of An- Therefore, from the measured All, one can get clear information of 
Ag/g with reliable calculation of dAapcF / dapcp- Q 
Furthermore, when Zi = Z2, we have another formula 

rfAa'^'^'^^ - rfAfx^^^^ - dAa""'""' ~ —Ag{r], Q^)dAapGF (H) 

9 

X {d{z, Q')D{z, Q^) + D{z, Q^)D{z, Q2) - 2 D{z, Q^)D{z, Q^)} . 

In this case, we can also define the double spin asymmetry All as follows; 

^ dAa^^^^ - dAa^-"^^ - dAa^'^' ^ ^9{v,Q^) dAapGF , . 

da^^^^ - da^-"-^ - da^--- g{v,Q^) ' dapGF ' 



The All defined in eg. (pUD and eq.(|T2|) results in the same physical quantity. From eq.(|T2p, we 
can also extract Ag/g as well as in eq . (|lO|) . 

3 Numerical calculations of the cross sections and the 
double spin asymmetry 

In this section, to see how the above formula works, we numerically calculate the double spin 
asymmetry for the large-pr pion pair production of pol-DIS. The spin-independent (spin- 
dependent) differential cross sections for producing hadrons hi and h2 are given byp 



d{A)a'^^>^^ _ d{A)a^,tF , c^(A)a^^^^j, ^^^^ 



dzid cos 9idz2d cos 92dxdy dzid cos 9idz2d cos 92dxdy dzid cos 9idz2d cos 92dxdy 



^Hcre we take account of the light quarks alone in eqs.(||)-(^). For large regions, heavy quarks might be 
generated for PGF process and QCD conipton process might also have contributions of heavy quarks. Even then, 
the ^LL is reduced to eq.® if D^* = D^* = D^" = . 



Each term in the right hand side of eq.([T3|) is written as 

dzid cos 9idz2d cos 92dxdy ' ' dzid cos Oidz^d cos 92dxgdy 

X E _ enDl^{z[,Q')D^^{z'„Q') + (1 ^ 2)} , (14) 

i=u,d,s,u,d,s 

-— ;y , , = E e A/,r^,Q2 — ^ . V , , 

dzid cos Uidz2d cos U2dxdy q-udsuds dZid COS UidZgd COS ti2dXqdy 

X {DjH^;,Q'K^(4Q') + (l-2)}, (15) 



where 



Tj = x + {I - x)tiT2 , Q = X y s , 

A = ( ^^^^ ) Zi , Z2= ( ^^^tZl j Z2 , (16) 



77x1X2 Scos^ cos^ i6'2 sin^ i(6'i + ^2 



with 



ri,2 = tan ^6^1,2 • 



Here we simply assume the scattering angle of outgoing hadrons 61^ 2 to be the same with the one 
of scattered partons in the virtual photon-nucleon cm. frame. This assumption might not be 
unreasonable if observed particles are light hadrons with high energy, s is the total squared energy 
of the lepton scattering off the nucleon. x, y and zi^ 2 in eqs.(p!3[)-(|T6D are familiar kinematic 
variables for semi-inclusive processes in DIS and are defined as 

_ _P-q _ P-Pi,2 

2F ■ q F ■ t F ■ q 

where g, F and Pi^2 are the momentum of the incident lepton, virtual photon, target nucleon 
and outgoing hadrons, respectively. The differential cross sections of hard scattering subprocesses 
with outgoing two partons having opposite in an azimuth angle for ^g i'qq and i Q ^ i' 9 Q 
at the leading order QCD are given as 

d{A)d-pGF{QCD) _ 1 a^as y{i] - x){l - r])'^ 



dzid cos 6id(f)idzi(^g-jd cos 92dxg(^q)dyd(f)£ 128tt^ {po ■ i)Q'^ x 

X Bie,,e2)eUlmM\laFiQCD). (17) 

with 

'{zi{l -r]) + {r]- x)} sin + {zi(g){l - v) + iv - a;)} sin 6*2 



^ = sin(^i + 62) 



sin 61 sin 62 



where Zi, zj and Zg are the momentum fraction of the outgoing parton i, i and g, respectively, to 
the incoming parton, and are given as||19|| 

- ^2 _ Ti 

n + T2 Ti + T2 

The amphtude |(A)M|p^p,(-Q^^-) in eq.([T7|) is given by ||2^ 

\i ■ Pi? + {i' ■ Plf + {i ■ P2? + {i' ■ P2? 



\M\Igf 



|M|^ 



16{i ■ £' 
128 



\QCD — 2 

for the spin-independent case and 



\^M\lcF 



\^M\lcD 



128 



{i ■ i') 



(Po ■ Pi)iPo ■ P2) 
(£-po)^ + (f •po)^ + (£-pi)^ + (f -Pi) 

{P0-P2){P1 ■P2) 



{i'-pi)'-{i-pif + {i'-p2f-{i-p2r 

{Po ■ Pi){po ■ P2) 

'{i-PoY-ii'-PoY-ii-PiY + ii'-Pi) 



(PO ■ P2)iPl ■ P2) 

for the spin-dependent case. The po, pi, P2 and i' are the momentun of the incoming parton, 
outgoing parton i, i, g and the outgoing lepton, respectively, and are written by 

Po = |P| (?7 , , , — rj) , pi = \pi\{l , sin 6'i cos 01 , sin 9i sin 0i , cos 9i) , 
P2 = b2|(l , sin6'2cos(0i - vr) , sin6'2sin(0i - vr) , cos 6^2) , 
i'^ = \£\{1 , sin6'f cos0£ , sin^^£sin0£ , cos6'£) , 



|P|(l-2x, 0,0,1), 



with 



\P\ = 
cos 9i 2 



\ 4x(l - x) 

^i,i(g)(l + ^-2x) -{rj-x) 



\pi,2\ = \P\{zij(^g){l ~v) + {v- x)} 

1 — xy 



P 



cosOe 



y 



1 — 2x + xy 
1 — xy 



in the cm. frame of the virtual photon-nucleon system|T9[. By using these formulas and 



newly analyzed pion fragmentation functions ||17||, we have calculated the spin-dependent and 
spin-independent cross sections of the large-pr pion pair prodcution and estimated the dou- 
ble spin asymmetry of eq.(|10]) at the energy of HERMES experiments. Here, we have taken the 
AAC[r[ and GS96p[ parametrizations at LO QCD as polarized parton distribution functions and 



GRV98|2T1 and MRST98|22| as unpolarized ones. At ./s = 7.25GeV, y = 0.75, > iGeV^ and 
W'^ > lOGeV^ with two different sets of kinematical values of 61 2 and Zi 2 for the produced pion 



pair, the calculated results of the spin-independent (spin-dependent) differential cross sections 
and All are shown as a function of 1] in Figs. 2 and 3, respectively. f\ From Fig. 3, one can 
see a big difference of the behavior of All depending on the models of Ag/g and hence, we can 
extract the behavior of Ag rather clearly from this analysis. 

4 Conclusion and discussion 

We proposed a new formula for extracting the polarized gluon distribution from the large-pr 
light hadron pair production in pol-DIS by making an appropriate combination of hadron pair 
productions. Since the double spin asymmetry All for this combination is directly proportinal 
to Ag/g, the measurement of this quantity is quite promising for getting rather clear information 
on the polarized gluon distribution in the nucleon. 

In this work, we calculated only the case of the large-p^ pion pair production. The same 
analysis can be applied also for the kaon or the proton pair productions by considering the reflec- 
tion symmetry along the V-spin axis, the isospin symmetry and charge conjugation invariance 
of the fragmentation functions as follows ITH 
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nK = 













Df = or = or = Df 



for the kaon pair production case and 



DP 






= Dl 


= D^ 


Jyp 




= Dl 


= Dl 


= Dl 






= Dl 


= Dl 


= Dl 













for the proton pair production case. 

^Here the polarized parton distributions were evolved from Qq = IGeV^ to any value. Though the initial Qg 
value of the GS96 model is taken to be 4GeV^ in original literature, we simply assumed their parton distributions 
to be in scaling for 1 < < 4GeV^. 



In general, it might not be so easy to precisely extract a physical quantity from a combination 
of data on the different kind of physical quantities. However, if we have enough data with high 
precision for wide kinematical region, it could be possible to determine the physical quantity 
concerned with rather precisely from those data, just as in the case of Bjorken sum rule in which 
two different quantities, Qi and (?", were combined. We believe that the formulas proposed here 
are suitable for those analyses, since we can expect to have rather many and precise data because 
of large cross sections for the light hadron pair production. The HERMES RICH detector at 
DESY which will provide good particle identification in the wide momentum range could allow 
to measure All for the processes discussed here. 

This work is supported by the Grant-in-Aid for Science Research, Ministry of Education, 
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p p 



(a) photon-gluon fusion (b) QCD compton 



Figure 1: Lowest order Feynman diagrams for the large-p^ hadron pair production. 




Figure 2: Combined spin-independent and spin-dependent differential cross sections defined at 
the denominator and numerator, respectively, of eq.(fO) as a function of rj at a/s = 7.25 GeV, 
y = 0.75 for the deep inelastic regions (Q^ > IGeV^ and > lOGeV^) with two different sets 
of kinematical values(^i^2, ^1,2) of the produced pion pair. 77 is the momentum fraction of the 
gluon and obtained from x by eq.(16). 




Figure 3: The rj dependence of All at ^/s — 7.25 GeV, y — 0.75 for two different sets of 
kinematical values(6'i^2, -^1,2) of the produced pion pair. Sohd hne and dotted hne are for AAC 
LO and GS96LO-C parametrization models, respectively. A.g/g itself is also presented for both 
parametrization models. 



